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I. INTRODUCTION 


A. . BACKGROUND 

The Captured Air Bubble surface effect ship presents 
a new concept in shipping for the surface force of the 
U.S. Navy. Within a few years these vessels, supported 
partially by a cushion of air and partially by the more 
conventional hull buoyancy forces, will be joining the 
Fleet. In preparation for that event it is of extreme 
importance for detailed studies to be undertaken so that 
the capabilities and limitations of such vessels may be 
fully understood. 

The term surface effect ship applies to a broad spec- 
trum of different types of vessels. Among them are the 
Air Cushion Vehicle (ACV), the Ground Effects Machine 
(GEM), the Semi-Submerged Ship (SSS), the Captured Air 
Bubble (CAB), and others. This study is concerned with 
the CAB surface effect ship. 

The Captured Air Bubble craft consists of two rigid 
Sidewalls which resemble catamaran hulls and flexible 
seals between the water surface and hull structure across 
the bow and stern of the craft. The enclosed volume is a 
plenum chamber into which air is blown by a lift fan system. 
The pressure in the plenum chamber acting over the area 
of the plenum provides a lift force sufficient to support 


a large proportion of the total craft weight. The remainder 
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of the supporting force is made up by conventional hydro- 
dynamic forces acting on the sidewalls, plus some lift 
forces due to the bow and stern seals and aerodynamic 
lift at higher operating speeds. 

In the “on bubble" mode of operation, at the higher 
operating speeds, approximately two-thirds of the total 
craft weight is supported by the air cushion, and therefore 
relatively shallow sidewall immersion is required to support 
the remaining weight. The shallow draft produces a lower 
hydrodynamic drag and thus higher speed for a given power 
than could be obtained for a displacement hull. 

The XR-3 is a CAB surface effect ship of approximately 
three tons displacement. The XR-3 was built by the Naval 
Ship Research and Development Center (NSRDC) in 1966. 
Following periods of testing by NSRDC, Annapolis Division 
NSRDC, and Aero jet-General Corporation, the XR-3 was de- 
livered to the Naval Postgraduate School (NPS), Monterey, 
California in March 1970 for the investigation of basic 
and advanced surface effect ship technology. 

In addition to the XR-3 there exist other CAB test 
craft. Among these are the Aerojet-General 100-A and the 
Bell Aerospace Systems 100-B, both of approximately 100 
tons displacement. These two vessels were constructed as 
prototypes of a 2,000 ton ship. To date speeds in the 80- 
knot range have been attained by these vessels. 

Under a U.S. Government contract, Oceanics Incorp- 


orated developed a digital computer simulation of the 
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CAB ship loads and motions for the Surface Effect Ship 
Project Office (SESPO) in Washington D.C. This loads and 
motions program with the mathematical model of the Bell 
Aerospace Systems 100-B with six-degrees-of-freedom has 
been in use at the W.R. Church Computer Center at NPS 
Since October 1972. 

In December 1973, LCDR. D.G. Leo, USN and LT. R. Boncal, 
USN produced a mathematical model with six-degrees-of- 
freedom of the XR-3 test craft derived from the above- 
mentioned simulation model of the 100-B (Ref. 1). The 
modular construction of the 100-B simulation program 
facilitated the re-modeling of the program where necessary 


to fit the design differences between the 100-B and the XR-3. 


B. OBJECTIVES 

It is the objective of this thesis to validate the 
non-linear, six-degrees-of-freedom, mathematical model of 
the XR-3 CAB test craft in a comprehensive manner by direct 
comparison between instrument-measured variables on the 
XR-3 and computed variables produced by the simulation 
programs. In the process of carrying out this objective, 
subroutine modifications were made to the simulation 
program to improve on the structure of the simulation model 
and justification is given in the sections that follow for 
the changes that were made. 

The data taken on the XR-3 measured in real-time, 


analog fashion using FM modulation was filtered, converted to 
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digital form, and processed in such a manner that it could 
be interfaced directly with the simulation program and a 
time history of true differences for each variable could 
be produced at each operating condition. 

The present study is limited to operation of the test 
craft in calm water conditions at the test site located 
at Lake San Antonio Reservoir, California. The range of 
operation of the test craft includes: 

1. Straight runs at various speeds with various 

locations of the center of gravity. The center 
of gravity was shifted by placement of two 125 

pound water tanks at different locations on the 
deck of the test craft. 

2. Straight runs at various speeds in an unladen 

condition. 

3. Runs made at various speeds with various rudder 


angles applied to validate turning performance. 
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II. COLLECTION OF DATA 


A. INSTRUMENTATION ON THE XR-3 

In order to support the continuing studies being under- 
taken at the Naval Postgraduate School on Captured Air 
Bubble vessels, the XR-3 has been instrumented with devices 
capable of measuring and recording time histories of many 
different variables for later, post-run analysis (Ref.2). 
The general configuration of the XR-3 is shown in Figure l. 
Table I shows a listing of the measurable parameters, the 
range of measurement, and the accuracy with which the 
Measurements may be made. The locations of the measuring 
instruments are shown in Table II. The measured data are 
recorded on a Pemco Model 120-B magnetic tape recorder 
which features 14 tracks for recording data and an edge 
track for voice annotation during measurement runs. Each 
of the 14 tracks on the tape recorder is capable of record- 


ing t 4 volt RMS with a tolerance of + 0.5%. 


B. POWER SUPPLIES 

The tape recorder is powered from a 110 volt 60 cycle 
auxiliary power unit, located on the aft weather deck as 
shown in Figure 1., through a Pemco a.c. to d.c. converter 
which provides a 26 volt d.c. power source. The three 
accelerometers for measuring vertical, lateral, and fore- 
and-aft accelerations are powered from a 12 volt storage 


battery through a d.c. to d.c. converter which provides a 
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Table I. 


Parameter 


Port Thrust 
Starboard Thrust 
Velocity 

Height 

Bowseal Press. 
Plenum Press. 
Sternseal Press. 
Pitch Angle 
Roll Angle 

Yaw Angle 

Pitch Rate 

Roll Rate 

Yaw Rate 

Rudder Position 


Longitudinal 
Acceleration 


Lateral 
Acceleration 


Vertical 
Acceleration 


*Height is used to compute Draft as 
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XR-3 Measurable Parameters 


Simulation Program Range 
Name 
THSTP 0-500 lbs 
THSTS 0-500 lbs 
VEL 0-40 kts 
None* - 2 feet 
PBARB 0-60 psf 
PBAR 0-60 psf 
PBARS 0-60 psf 
THETAR £15 degs 
DPHI 715 degs 
DPSI * 180 dees 
PDEG 730 degs/sec 
QDEC +30 degs/sec 
RDEG “30 degs/sec 
DELRS * 45 degs 
None - 0.2g units 
ACCLAT - 0.2g units 
ACCEL (3) 


Accuracy 


I+ 


5 1D 


I+ 


510 


i+ 


1 kt 


I+ 


O.lin 


t+ 


O25 pst 


I+ 


0.5 psf 


I+ 


0.5 pst 


+ 


0.5 deg 


1+ 


0.5 deg 


1+ 


0.5 deg 


20.5 deg/sec 
te 2. y 
-0.5 deg/sec 


t0.5 deg/sec 


not available 
not available 


not available 


- 0.8 to +1.2g not available 


explained in text. 
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Table II. Location of Sensors 


Sensor Distance Forward Distance from 
of Transom Center Line 
Pressure 
Transducers 
Bow Seal 18'4" . € 
Plenum 12'yan € 
Stern Seal 1G" 
Velocity otto" € 
4°93" down from 
deck 
Height ena" € 
6" down from 
deck 
Gyros 
Pitch 1276" 2°*4s" Port 
Roll 138 2°*4s" Port 
Yaw to" 2743" Port 
Port Thrust @) 5'6" Port 
Starboard Thrust 0) 5'6" Starboard 
Vertical 
Accelerometer 13°63" 1°95" Port 
Lateral 
Accelerometer 13 G" 179" “Port 
Longitudinal 
Accelerometer 13°53" 1°83" Port 
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28 volt d.c. supply. The three gyros which measure pitch, 


roll, yaw, pitch rate, roll rate and yaw rate are powered 


by the 12 volt storage battery. The gyros also receive a 


1 volt d.c. excitation signal from the 110 volt a.c. power 


supply fed through an a.c. to d.c. converter. 


The height sensor mounted forward of the bow is 


powered directly from the 110 volt a.c. auxiliary power 


unit. The height measurement signal is converted to a 


71.0 volt d.c. Signal by a height sensor conditioner. The 


velocity probe, 
sensor, extends 
testcraft. MThe 
probe is the 12 


converter. The 


mounted on the same assembly as the height 
into the undisturbed water forward of the 
28 volt d.c. power source for the velocity 
volt battery through the d.c. to d.c. 


rudder position measurement is made by means 


of a load cell which receives a 12 volt d.c. supply from 


the battery. 


The two load cells for measuring port and starboard 


thrust and the three pressure transducers for measuring 


bow seal, plenum, and stern seal pressures are all powered 


by a 28 volt d.c. supply from the 12 volt d.c. storage 


battery through 


the d.c. to d.c. converter. All five 


units also receive a 5 volt d.c. excitation signal from 


an amplifier-transducer package. 
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TIT. DATA PROCESSING 


A. ANALOG TO DIGITAL CONVERSION 

To provide measurements for validation of the mathe- 
Matical model, data-gathering runs were made with the XR-3 
testcraft on four different days at the Lake San Antonio, 
California test site. The data was recorded on magnetic 
tape using FM modulation, then processed into a suitable 
form for comparison with corresponding values computed in 
the XR-3 simulation program on the IBM 360/67 Digital 
Computer. 

The data was converted from analog to digital form 
on an XDS-9300/CI-5000 hybrid computer at the NPS 
Computer Laboratory. From previous experience in data 
reduction by the XR-3 Research Group in the Aeronautical 
Engineering Department at NPS, it was found that filtering 
of the analog signal was required to remove unwanted high 
frequency noise from the recorded data. To provide this 
required filtering six Krohn-Hite digitally-tuned filters 
were assembled so that six tracks of recorded data could 


be processed simultaneously. 


B. DESCRIPTION OF EQUIPMENT 

The analog tape containing raw data from XR-3 test 
runs was mounted on a Honeywell Tape Recorder/Reproducer, 
Model 96 for playback. Connections were made through 


trunk lines to lead the output from the six selected tracks 
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to the COMCOR CI-5000 Analog Computer as shown in Figure 2. 
Each channel of recorded data was then passed through an 
amplifier with a gain of 5 before being led out of the 
CI-5000 to one of the filters. The filters were operated 
in the low pass-maximally flat mode which provided a four- 
pole Butterworth response with an attenuation of 24db per 
octave from the cut off frequency. Cut off frequencies 
were set at 10 Hertz for all channels except port and 
Starboard thrust which were set at 1.0 Hertz. Each 

Signal at the output of its respective filter was led back 
into the CI-5000 for another amplification by a gain factor 
of ten. A Brush 8-Channel Recorder was also patched into 
the system at this point so that each of the six selected 
data tracks could be continuously monitored. Figure 2 
shows the signal flow patch through the hybrid computer 
system for analog-to-digital conversion. 

For analog-to-digital conversion, digital-to-analog 
conversion, and other applications, the CI-5000 Analog 
Computer and the XDS-9300 Digital Computer are intercon- 
nected through trunks. For A-to-D conversion in this 
application a sampling frequency of 200 Hz was chosen. In 
steady-state calm-water runs with the XR-3 none of the 
measured parameters can be expected to vary at frequencies 
greater than 10 Hz. Therefore, selection of a 200 Hz 
sampling frequency more than satisfies the requirements 


of the Sampling Theorem. 
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The first step in the A-to-D process was to zero all 
possible bias errors in the system. Calibration steps 
are recorded on separate channel strips on the analog 
data tape which indicate zero and maximum signal readings 
for each sensor. Bias potentiometers at the input of the 
first stage amplifier on the CI-5000 were adjusted to 
properly zero all channels during the pre-test run cali- 
bration of the XR-3 sensors. Each signal was monitored 
on a digital voltmeter during the calibration procedure. 

Execution of the FORTRAN program listed in Appendix A 
on the hybrid computer system enabled selective conversion 
from analog to digital form of specific time periods of 
analog data. Continuous monitoring of the voice edge track 
on the data tave provided a reference as to what the XR-3 
conditions were at all times. After each file of data had 
been digitized and written on the digital magnetic tape, 
the tape was re-wound and that file most recently recorded 
was run through the digital-to-analog conversion. The 
D-to-A signal was monitored on the Brush Recorder and com- 
pared with the analog signal to verify the validity of the 
conversion process. The finished product of the A-to-D 
conversion was a reel of 7-track digital tape consisting 
of many files of digital data. Each file on the tape 
contained a time history of six selected variables measured 


for a specific operating condition on the XR-3. 
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C. INTERFACE WITH IBM SYSTEM 360/67 

Introduction of the 7-track tapes into IBM System 
360/67 at the W.R. Church Computer Center at NPS entailed 
some further processing of data. In the XDS 9300 a binary 
number is represented by considering its octal representa- 
tion, made up of 3-bit digits, whereas in the IBM 360 
the same number is represented by hexadecimal, made up 
of 4-bit digits. As an example of this difference, a 
decimal number (668), 5 is (1234), octal and (290) 4¢ 
hexadecimal. It was necessary, then, to write a FORTRAN 
program to read each record in a file from the 7-track 
tape, make the conversion from octal to hexadecimal, store 
the data in hexadecimal form on a temporary storage disk, 
and to print out a portion of each record in decimal form. 
A listing of the program is included as Appendix B. The 
program uses the pre-compiled assembler language subroutine 
FORM to convert the data to its hexadecimal representation. 
Figure 3 shows the flow of information through the IBM 
System 360/67 computer in this validation. 

The printed output from the conversion program is re- 
presentative of the data in its original measurement form 
multiplied by a gain factor of 50 in the A-to-D conversion 
process. At this point in the process the numbers printed 
out were compared with the signals recorded on the Brush 
Recorder during A-to-D conversion to insure that the 
accuracy of the data had been preserved. This was indeed 


the case with the exception of some small amplitude, 
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high-frequency fluctuations in the digitally reproduced 
data which did not appear to be present on the Brush 
Recorder tapes. 

The final phase in data processing was to recover the 
Measured data into the units which correspond with the 
measurement, i.e. pounds of thrust, degrees of pitch angle, 
knots of velocity, etc. to facilitate direct comparison 
in the XR-3 simulation model. To accomplish the recovery 
of data it was necessary to read the converted data from 
the temporary storage disk, multiply each set of variables 
by the inverse of the appropriate scale factor used when 
recording the data on the testcraft, divide out the gain 
factor of 50 introduced in the A-to-D process, and store 
the result on a temporary disk for later access by the 
XR-3 simulation program. 

The height measurement made on the XR-3 during data 
collection runs does not have an immediately available 
counterpart in the XR-3 program as indicated in Table 1. 
A very close approximation to draft can be obtained from 
the height measurement through suitable scaling, however. 
The calculation of draft in the XR-3 simulation program 
.is made with reference to draft at the center of gravity 
of the craft. The height measurement is made at a sensor 
located 191" forward of the reference center of gravity 
and 6.375" down from the deck level of the test craft 
(see Figure 1). The measurement of the height of the 


sensor above the surface of the water will vary as the 
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craft pitches up and down at the bow. Since pitch angles 
are, for the normal operating range, less than five degrees, 
the small angle approximation may be made for the sine of 
pitch angle. The draft measurement is synthesized by 
solution of the following relationship: 

DRAFT = SH + Z - HT + D*SING 


where DRAFT = Water height above keel at the center of 
gravity 


7.145" 
Z = Distance of CG above keel = 30.48" 


SH = Height of sensor above CG 


HT = Measured height of sensor above water surface 
D = Distance of height sensor forward of CG = 191" 
D*SIN@ = Variation of measured height due to pitch angle, 
A FORTRAN program was written to accomplish the process 
of re-scaling to measured units and to provide the ciose 
approximation to a draft measurement. A time index was also 
included in the program to facilitate plotting each variable 
versus time on a CALCOMP plotter. This feature proved to be 
valuable as a final step in the process because each variable 
could be closely examined in graphical form prior to acces- 
sing the data by the XR-3 simulation program. It was at 
this point in the procedure that some difficulties were 
encountered in early attempts at signal processing, as 


described in the following section. 


D. FREQUENCY COMPONENTS ENCOUNTERED IN DATA 
Examination of the CALCOMP plots of the digitally 


reproduced measurement data revealed significant small 
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amplitude, high frequency noise impressed on the measure- 
ment signal. The noisy signal was not immediately evident 
on the Brush Recorder tape, however on close inspection it 
was apparent that some noise was passing through the filters. 
Figures 4,5,6, and 7 are representative CALCOMP plots of 
processed measurement data of thrust, pitch angle, velocity 
and plenum pressure respectively. These graphs are in- 
cluded to show the effect of noise on the processed data. 
The source of noise present in the measurement data 
has not been established and can only be rationalized to 
be a combination of several possible environmental factors. 
First of these is vibration aboard the XR-3 testcraft from 
the twin engine propulsion system and fan drive system. 
Second, the testcraft is rarely run in completely glass- 
calm water conditions at the test site. Surface disturbances 
are nearly always present due to numerous pleasure craft, 
wind-induced ripples on the lake, and occasionally, the 
wake of the testcraft itself. Such disturbances can have 
a pronounced effect on the motion of the XR-3. Finally, 
60-cycle power supply interference can be assumed to be 
present due to the antenna effect of unshielded wiring 


both in the testcraft and in the NPS Computer Laboratory. 


1. Elimination of Unwanted High Frequencies 
Recognizing that the noisy data rendered it all but 


unusable for any meaningful attempt at validation, methods 


were sought to eliminate the noise and to obtain a "clean" 


29 


mee" « 


i : 
- 0 


tnayecl saoto mo Sevewo veqat tobr09 8h dew 


Wy i f 
ouit aolegeqg asw sclom amos tae davseqin 

TMOOIAD ov itattarzerde: sta ) bre ita a 
Aas foviq ,tauaiy Io stab tasnotieasm hea ae . 


16 eiqaety saont .yievitosqaer esswess7q mune 
- ip 


; ' 
_% ar " ~ 
ie | “ fe . ¥~Gil£o ol J ans or 
i Og : Vmee2sitd « vit 20 SO 1N0G 
’ 
+ | he > $ 
od Yiile f berulidates head oe 
a 
j it . 
i 7 i. Eke | Se | 4 if 
* j {4 \ +, ; } os 
i P - =e 
f - ‘, fsiee e? 
ett 1071q otto 
J Lt as LoT SS 3 
* 
5 1 
rf CLS LDCS 
‘ : t : pT Th 


boa enel alt mo. es 1a it beouba f- 
-tioetl Ptatotaet aft ‘to 


b o4R -oitom att 00 teste bean) 


stistil ylqque aawog @ 


iw fx 0 osetia sonnets op of seb tas 


nese fet gan ott ti bos dtevoraed ne 


inthe date pstagwy, Xo sobiariih ie. « 
ion ated “yaabe ett tart? artshignoon 


eS ie. met 
iLisy va rymate ee: ‘sot 


a se 


cr 


Dan 


Ove 


005 


Oot 


QU 


Thrust 
(100 
lbs/inch) 


NZ 


Ooi 


Figure 4. 


Gel tg 


Patent a A a eed ey oy AA HAN nL Ay I TN RY QU NE 


= mF at 
vu2 You 945 


Time 
(2 sec/inch) 


Plot of Thrust Vs. Time 


30 


- 


Ag J 
nit ey, peurtt to | off, 
d 7 a 4 


Fe 9 
ayia 


uae | 


- 
is 
: Ty) 


Pitch 
Angle 


(0.50 


inch) 


¥. 
a 


1” 


yw 
a 


019 


Ws 


~ 
. 


we 


_s 


Figure 5. 


wif! | Hil 

ete gett a iy 
Hh fe a aa in Hit 
yet a1! “Ti Hy st sia 

‘| Nip ] i i 

ike t ua! 

Mi 
Ke O24 acs 


Plot of Pitch Angle Vs. 


ie fi 


Time 


th, 


hi Hingis 


vty, 
sila! : i 


‘agai 
GCE rei | 
Time 


(2 sec/inch) 


035 


030 


920 


Velocity 
(5 kts/ 


inch) « 


A IAA OY RA aA EA A RAEN AAS MN RARE 


low) 


Sac Uc2 Qu; 055 oo8 20° 


Time 
(2 sec/inch) 


Figure 6. Plot of Velocity Vs. Time 


32 


— } 


Cae 


dbs xy Nr tootey 


939 


Ly NUUOUUUOLUPE auaidll Malacca cen ‘ LRA Ag uaa A CAMORPARA yah 
i Re aia ail ws "even 


oO 
N 
ao 


Plenum 
Pressure 


(5 lbs/ 
inch) .. 


‘ 
Led 


019 


oct 


Time 
(2 sec/inch) 
Figure 7. Plot of Plenum Pressure Vs. Time 


33 


=. . ao 
Sec vce QGi 965 RRL} S 


signal and, at the same time, to preserve the validity of 
the actual measurement. The first effort was to filter 
at a lower frequency. The result was displeasing for, 
while the magnitude of the impressed noise was somewhat 
reduced, lower frequency oscillations which represented 
actual craft motions were nearly eliminated. The next 
possible solution was to employ two filters in cascade, 
thereby obtaining a low-pass filter with 48 db per octave 
attenuation above the cut off frequency. This procedure 
accomplished the main objective of removing virtually all 
the noise and passing only signals representative of the 
craft motions. The disadvantage of this method was that 
only three instead of six measurements were now possible 
because of the hardware limitations and synchronization 
became an insurmountable problem. 

The third possibility of removing the noise from 
the measurement signal was to proceed on the assumption 
that the noise was gaussian and white and to employ a 
smoothing or time-averaging technique. The technique em- 
ployed was to process the data as originally set forth and 
then, in the final scaling process on the digital computer, 
to obtain a moving time average. The algorithm implemented 
on the digital computer was to calculate the sum of the 
first ten points in the time history of a measured para- 
meter, divide the sum by ten and call the resultant figure 
the first point in the time history. Next, a counter was 


incremented by five and the process was repeated to obtain 
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the next point. In this manner the data is reduced from 
200 samples per second to 20 samples per second. The over- 
lapping smoothing technique tends to preserve low frequencies. 
By comparison there were insignificant differences (less 
than two per cent) between this smoothing process and the 
aforementioned cascade filtering. The advantage realized 
was a more economical use of available equipment (filters) 
and a great savings in time required for data processing. 

A listing of the program used for re-scaling the 
data to its original measurement form and including the 
smoothing routine is included as Appendix C. Figures 8, 

9, 10 and 11 are CALCOMP plots of smoothed measurement data 
of thrust, pitch angle, velocity and plenum pressure 
respectively. By comparison with the measurement data 
plotted in Figures 4,5,6, and 7, it is evident that the 
moving time average smoothing techniaue has preserved the 
measurement information of value and has eliminated the 


impressed noise. 


‘Tes = 
ee | ho 
‘Sr? 


oe 9 eran t wel 
19Go| Of Rian suplestont ge Ltd oe sat 


> 


» tassttingtaal atew eayity roa uscfi 
= 


itt Game elett meewi ed (¢030 aaq) ah 


> “ “ 
0 } c re obs 
- ~eis : 
d t Be he ow 
‘ 
a 
bat + %o 
=J09R & 
i « 7 
; soe a 


rtova oo , 


| sora! tne 


09? 


036 


Out 


Thrust @ 
(100 lbs 
inch) 
3 
Figure 8. 


Plot of Thrust Vs. Time After 


36 


ues aes] 


Time 
(2 sec/inch) 


Smoothing 


D| 


jx-~atsth ith get th onl siete ACC ety st aly 


jee a ; eee: in - 
Ths emit 6 
font \oeu.X}- 
® va z i 
anids Cage 19th ont? ey Sane 
: ee 7 


; a aud a. i 


. u —_ 7 e 


Pitch 
Angle 


(0.5° 


incn) « 


O10 


re) 
C3 
2 


Ww 


=a 


rs An7e ace 


Sos aG2 <4 955 oF oid 


Time 
(2 sec/inch) 


Figure 9. Plot of Pitch Angle Vs. Time After Smoothing 


37 


} 
a i 
ye 5 


oe Tete Lae a 


; a 
a Dany cy : ’ 
i. - 


a 


re: 
ae 


020 


g 
Velocity 
(5 kts/ 
. w 
inch) da 
a 
=4 
oO 
re) 
$ 
qQ 
pay 
ate a2 033 025 Soy aig 
Time 


(2 sec/inch) 
Figure 10. Plot of Velocity Vs. Time After Smoothing 


38 


ee ae 


‘tian \oon Sp 7. ue 

vid cade tatiana vite 
. (Oh «tore 

‘ i ; ee id 

-_ ice ai | 

Poe LL tow 


039 


Plenum 
Press. 


(5 10s 
inch) f 


019 


vob 


10 


vee 


Figure 11. 


=~ - 
weve vues 


Plot of Plenum Pressure 


29 


V 


Oe 


Time 


aa) i) 
v3 wou 


Time 
(2 sec/inch) 


After Smoothing 


7 
: 
7 
J 
; 


IV. INTERFACE WITH XR-3 SIMULATION PROGRAM 


A. LOGIC ADDITIONS TO SIMULATION PROGRAM 

The XR-3 Loads and Motions simulation program was 
modified to the extent that, on execution of the program, 
the processed measurement data could be accessed if desired 
and on-line computation of differences could be made. ‘The 
output from the program would then be a time history of 
the differences between computed and measured parameters. 
This entailed the addition of another option switch in the 
input data which is | by subroutine INCON. If the option 
Switch is on or equal to one, subroutine INCON reads the 
measurement data into core from the disk where the data 
is stored. If the option switch is off or equal to zero, 
no reading of measurement data occurs. 

It was decided to restrict validation or data compari- 
son runs to ten seconds of simulation time because of core 
storage and computer time allocations. This decision set 
the size of the measurement data arrays at 200 elements. 

A COMMON block was added in subroutine INCON and RHS so 
that the measurement data could be available in subroutine 
RHS. As with other COMMON blocks, the feature to zero all 
elements of the added COMMON block was included in the 
Main program. The following XR-3 test craft variables 

and the associated simulation names were selected for 


validation: 
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Total Thrust THST 


Velocity VEL 
Height DRFT 
Pitch Angle THETAR 
Roll Angle DPHI 
Plenum Pressure PBAR 
Yaw Angle DEPSI 
Rudder Angle DELRS 
Vertical Acceleration ACCEL(3) 
Yaw Rate RDEG 


The necessary logic was added in subroutine RHS such 
that the appropriate element with respect to time in the 
measurement array could be compared with the computed 
variable and a difference formed. The resultant differ- 
ences were written into a temporary file for output by 
subroutine COLFIL after the simulation run reached the 
specified final time of ten seconds. Computation of 
differences was arranged such that positive values indi- 
cated that the computed parameter was greater than the 
measured parameter for any given sample time. The print 


interval for output purposes was set at one point every 


0.05 seconds to correspond with the measurement data rate. 


B. METHOD OF THRUST INPUT AND OBSERVED DIFFERENCES 


The first step in the validation process was to 


initialize the mathematical model at a particular operating 


point with the expected trim conditions of velocity, pitch 
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angle, draft, and plenum pressure. The feature exists in 
the simulation program to provide a thrust map of up to 

25 points. Additional logic was included to select every 
fourth point in the 200 element array of measurement 
thrust to be applied to the mathematical model as an input 
variable and a ten second simulation proceeded to produce 
computed velocity. All other measurement parameters except 
velocity were ignored during initialization runs. The 
difference between computed velocity and measured velocity 
was taken as an output from the model. Following initial- 
ization of the model, another computer simulation was made 
using computed thrust as necessary to hold velocity con- 
stant. All measurement data for the run was read in and 
differences between computed and measured parameters were 
taken as output from the computer model. 

For validation studies using measurement data with 
very nearly constant thrust applied on the testcraft, the 
foregoing procedure could be shortened by including only 
the last phase. Provided that the initial trim conditions 
were known for input to the computer model, no loss of 


accuracy resulted in the process of validation. 


C. METHOD OF RUDDER AND THRUST INPUT TO THE MODEL 
Similar to the thrust map provision is a provision 

for a rudder map input to the computer model. Up to 25 

points may be included in a rudder map and the necessary 


logic was included in the program to select every fourth 


42 


‘, 7 
‘bohulowts Baws 


a Lelion Leo ti aaentet aid 


e80N'Ig ivivsiuais bnosem net 


Temsiag Ihe Jea29m reo LEA? wt 


cit aobysrilalvini atiavh begorgs #19W Wee 
a 
? i] 


beivaeen bie YFRooley betvques neawi ad’ eon 


: tivo {ie% , Letom mort fvugqive Ne ee 


} fz Rica 

10 £3 to atsdouqg/ edt GF bed Luee? ¥ 
i hd 2) n 
UWNr Teun Ona. ROR . a0! 

ad rah 
f nels i 2 is (Zi iO” 8 B31) ravads ate 08 | 
ae ‘WJ Dog a i ach ot P > 
A al ae ' a wy 


point in the 200 element array of measured rudder angle 

for input to the mathematical model. The objective here 
was to apply the same inputs to the mathematical model as 
were used on the testcraft along with the initial trim 
conditions for each case. This method proved to be of 
little real value because a steering cable failure to the 
port propulsion engine on the testcraft rendered the rudder 
position measurement unreliable. To substitute for direct 
steering control on the port engine, the starboard engine 
was mechanically linked to the port engine by means of a 
long, rigid bar. Thus, the starboard engine controls the 
port engine in the rudder position and alignment of thrust 
vectors on the testcraft. Virtually all precision in rudder 
angle measurement is lost due to unavoidable play in the 
linkage and lack of direct steering control over both 
engines. 

Validation studies involving turning maneuvers with 
the XR-3 testcraft and the mathematical model were carried 
out by matching yaw angles and yaw rates for a spectrum of 
operating conditions and by noting the differences between 
computed and measured rudder angles for each operating 


condition. 
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V. PROGRAM SUBROUTINE MODIFICATIONS 


A. REPLACEMENT OF CENTER OF PRESSURE CURVE 

Newman and Poole reported in Ref. 3 the existence of 
a pressure wave and the attendant gradual deformation of 
the surface of the water due to a moving pressure distribu- 
tion with regard to CAB type hulls. NSRDC found that, in 
towing tank tests of CAB scale models, the craft tends to 
pivot in pitch motion about a point approximately 75-80% 
of the wetted length of the sidewall forward of the tran- 
som (Ref. 4). Based on these findings, it has been 
postulated that the center of volume within the plenum 
chamber moves aft with increasing speed. 

Accordingly, the caiculations made in the XR-3 sim- 
ulation program to account for a shifting center of pressure 
were removed from subroutines RHS and WAVES. Computations 
were inserted in the program to account for the addition 
of an incremental wedge of volume to plenum volume which 
varies with craft speed. What follows is a description 
by subroutine of changes which were entered and the basis 
for the inclusion of the changes. 

1. Subroutine INCON 

A previously unused space on the input data card 
number 00701 was used to read in the location of the pres- 
sure wave pivot point, measured in feet and tenths of 


feet forward of the transom. The statement XPWV = TEMP(3) 
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was inserted in Block 07 in the INCON subroutine. A value 
of 17.2 ft. was used for XPWV in this study. 

The relationship for wave drag produced by the move- 
ment of the pressure bubble over the water surface is, 


from Ref. 5: 
Aw 
F —_ 
pel 
The same formula is mechanized in the simulation program as: 


— : : 
bubpie = = (P,, PS) C 


FXPWAV = -PBBAR * CF * PWVCON 
where CF = 0.37 / ((U/FNCON) **1.5655981) and PWVCON = 
4,0* WEIGHT /(RHO *G* XLBW). The coefficient 0.37 and the 
exponent 1.5655981 are correct for the length to beam ratio 
of the XR-3 as determined from Ref. 3. The variable name 
WATSLP is given to the slope of the water surface inside the 
plenum and is calculated as: 

WATSLP = PBBAR * CF * PWVCON / WEIGHT. 
WATSLP represents the angle at the apex of the wedge-shaped 
increment of volume to be added to plenum volume. Since 
the ratio of PBBAR*CF*PWVCON ( i.e. bubble drag) to WEIGHT 
will always result in a small number for speeds above hump 
speed, the small angle approximation may be used without 
loss of accuracy. Therefore, the volume increment due to 
the water slope to be added to plenum volume is given by 

0.5 * WATSLP * XL * AB 

where XL is the length of the plenum chamber and AB is the 
area of the plenum at the water surface. The quantity 0.5 
* WATSLP * XL is then the vertical cross sectional area 


of the wedge-shaped volume increment. 
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Two other calculations were added in subroutine 
INCON for use in other subroutines during the simulation. 
The distance from XPWV from the forward extremity of the 
plenum is calculated as 
XLXPWV = XLBW - XPWV 
where XLBW represents the plenum length at the water 
surface. Also the distance between the craft center of 
gravity and XPWV is calculated as 
XPWVXS = XPWV - XS 
where XS is the longitudinal center of gravity measured in 
feet forward of the transom. 
2. Subroutine BOWSL 
In the computation for ELSKI, the wetted lenght of 
the bow seal, an incremental length is added to the wetted 
length of the seal. The increase is due to the water slope 
and is calculated simply as 
XLXPWV * WATSLP 
and added to the existing expression for ELSKI. 
3. Subroutine STNSL 
Similar to the wetted length correction in BOWSL, 
an incremental correction is computed and applied in 
subroutine STNSL. Here the correction is opposite in 
Sign and is applied to the wetted length of the stern seal, 
ELSKI as 
-XPWV * WATSLP. 
4, Subroutine RHS 
The water slope correction is computed in subroutine 


RHS as a continuous variable over the period of any given 
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Simulation time. Recognizing that WATSLP is calculated as 
the ratio of bubble wave drag to total craft weight and 
that bubble wave drag varies with velocity and plenum 
pressure, WATSLP must be considered as a continuous 
variable. Accordingly, the following computation was in- 
troduced in RHS: 

WATSLP = -FXPWAV / WEIGHT 
The plenum volume is corrected for the water slope on each 
pass through subroutine RHS by the addition of a volume 
increment as in subroutine INCON: 

+ .5 * WATSLP * XL * AB. 

5. Subroutine SIDEWL 
The addition of a wedge-shaped volume to plenum 

volume required an alteration to the computation of the 
drafts along the sidewalls inside the plenum chamber. 
Displacement of water is calculated as 

PBHEAD = PBAR / (RHO * G) 
where PBAR is the instantaneous plenum pressure. The im- 
mersion depth on the inside of each sidewall section is 
calculated in the relation 

DDIN = DD - WATSLP * (XPWVXS - XX(,K)) 
where DD represents the uncorrected water level height 
above the keel along each sidewall section and XX(J,K) 
is the distance from the transom of the center of each of 
the eleven stations along the sidewall. 

A revision of the method of computing sidewall gaps 


at any given station was also necessary. A logical IF 
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statement was included to test each value of DDIN after 
it was computed. If DDIN is less than zero another logical 
IF statement must be tested to determine if a gap exists 
between the bottom of the sidewall (keel) and the water 
surface. If DSW (J,K), or the actual sidewall immersion 
depth, is less than PBHEAD, then a gap for that sidewall 
section is computed as 

GAP (J,K) = - DDIN * (1.-(DSW(J,K))/ PBHEAD) 
where the computed GAP for any given station is in units 
of feet. 

Further background information of use in this 


study is contained in References 6-11. 


B. THRUST VECTORS IN SUBROUTINE PROP 

The propuision system in the XR-3 consists of two 
forty horsepower outboard motors. The port screw has a 
right hand rotation and the starboard screw has a left 
hand rotation. Steering control is obtained by rotation of 
the motors in the horizontal plane. fMThus, there is some 
rudder effect from the foil-shaped lower motor housings, 
however the predominant turning moment results from re- 
alignment of the propulsion thrust vectors when rudder angle 
is applied. In addition to the main thrust vector, directed 
in the line of craft motion for zero rudder angle, there 
exists a side thrust vector for each propeller. Since the 
propellers are counter-rotating, and due to the direction 


of rotation of each propeller, the side thrust vectors are 
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oppositely directed and, for zero rudder angle, are 
directed inboard. 

The assumption is made that both port and starboard 
propulsion motors rotate together through the same rudder 
angle. It is also assumed that the magnitude of the main 
and side thrust vectors are independent of rudder angle. 

The propeller thrust diagram shown in Figure 12 was 
used to verify and correct the force and moment equations 
for the XR-3 propulsion system. The figure displays the 
Main thrust and side thrust vectors for both propellers 
positioned with right (positive) rudder angle§, positive 
pitch angle6, and negative (port side down) roll angleg. 
The thrust vectors act through the port and starboard pro- 
peller positions, XP, YP, and ZP from the craft center of 
gravity in the Body Reference coordinate system. 

The forces associated with the propeller alignment of 
Figure 1 are shown vectorially in Figure 13 in the Local 
Level coordinate system after having been corrected for 
pitch and roll angles. The assumption is made that pitch 
angle and roll angle are limited to small angles such that 
siné=@ and cos@= 1, sing =¢ and cos# = 1. The follow- 


ing force equations result: 


T1X = THST1 * cos& 
T2X = THST2 * cos§ 
ST1X = - STHST1 * sind 


ST2X = STHST2 * sin$ 


T1Y = THST1 * sin $ 
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T2Y = THST2 * sin& 

ST1Y = -STHST1 * cos § 

ST2Y = STHST2 * cos § 

T1Z = THST1 * cos§ * ® 

T2Z = THST2 * cos& *® 

ST1Z = STHST1 * cos * 

ST2Z = STHST2 * cos§ *9 

Summing the forces in each of the three corrdinate 
directions yields the following result: 


FXS = T1iX - STIX 


FXP T2X + ST2x 

FYS = - T1Y - ST1Y 

EYP =-T2Y++ SiT2y 

FZS = - 712 + ST1Z 

FZP = -12Z - S12zZ 

The total force equations are found by summation of 
the force equations for the port and starboard propulsion 
systems in each of the three coordinate directions as 


follows: 
FX 


FAP + FXS 

FY = FYP + FYS 

FZ = FZP + FZS 
Corrections were made to the force calculations of sub- 
routine PROP from Reference 1 on the basis of the foregoing 


discussion. 
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VI. DISCUSSION AND EVALUATION OF RESULTS 


A. OBTAINING AGREEMENT AMONG VARIABLES 

Computer simulation runs in the initial stages of this 
study revealed that plenum pressure was consistently ap- 
proximately 2.5 psf higher than that recorded on the XR-3 
testcraft. Among other variables, thrust on the testcraft 
was higher than the computer model indicated for any given 
speed and calculated draft was greater by approximately 
three inches than the computer model produced. The decision 
was made to first obtain agreement in plenum pressure 
between the mathematical model and the testcraft. 

The bow seal leakage area input to the mathematical 
model was listed in Ref. i as one of the parameters not 
known to be exact. The estimate of leakage area was re- 
vised trom 0,00: 6q,. £t. to O0V10 sq. ft. ‘This slight alter- 
ation had the desired effect of reducing plenum pressure 
to agree with that measured on the XR-3 testcraft. An 
unexpected but welcome additional effect was nearly perfect 
agreement in thrust and draft. Aside from the alterations 
to the structure of the mathematical model there were no 
further changes made to the model during the course of 
the present study. 

Ty 


Straight Runs Made with Different C.G. Positions 


Validation runs were made interfacing data measured 


on the XR-3 testcraft with the computer model. The data 
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presented in Table III is a summary of average differences 
calculated between time histories of actual craft motions 
and simulated craft motions for the thrust, velocity, 

pitch angle, plenum pressure and draft variables. Also 
included is percentage differences for ease in qualitative 
comparison. It should be noted that percentage differences 
in pitch angle measurement appear to be quite large, how- 
ever in most cases the difference is within the measurement 
accuracy limitation of = 0.5 degree aboard the testcraft. 

The total craft weight on the measurement runs 
listed in Table III was 6270 lbs. In addition to the pilot 
there was a passenger in the co-pilot position and 250 lbs. 
of water in plastic containers aboard used for changing 
the location of the center of gravity of the XR-3. 

The data presented in Table IV is a summary of 
average time history differences with a fixed center of 
gravity on the XR-3. On the measurement runs listed in 
Table IV the total craft weight was 6050 lbs. Vertical 
acceleration was substituted for draft in this series of 
validation runs. 

Figures 14, 15, 16, 17, and 18 are typical of the 
graphical output from the XR-3 simulation program arranged 
to display the time histories of differences between com- 
puted and measured variables for this study. 

2. Runs Made with Turns 

A series of eleven validation runs were made with 


various rudder angles at various speeds. A summary of the 
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Figure 14%. Plot of Thrust Difference Vs. Time 
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Figure 17. Plot of Plenum Pressure Difference Vs. Time 
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Figure 18. Plot of Draft Difference Vs. Time 
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average differences calculated between actual craft motions 
and simulated craft moticns during the turning maneuvers 
for thrust, velocity, pitch angle, roll angle, yaw angle, 
yaw rate and rudder position variables is presented in 
Table V. 

The relatively large differences in rudder angles 
in many runs is attributed to the lack of confidence in 
an accurate rudder position measurement. In the future 
when direct steering control of both propulsion motors is 
restored, it is recommended that confirmation be made of 
the results presented here with regard to rudder position 
and other variables in turning maneuvers. 

Figures 19, 20, 21 and 22 are typical of the graph- 
ical output from the XR-3 simulation program which display 
the time histories of differences in roll angle, yaw angle, 


yaw rate, and rudder angle, respectively. 


B. LIMITED TO CALM WATER CASES 

It is recognized that much of the measurement data 
obtained at the XR-3 test site at Lake San Antonio, Cali- 
fornia contains small perturbations due to surface 
disturbances on the lake. No attempt has been made to 
measure the surface disturbances or their effect on the 
testcraft. Validation of the mathematical model in the 
presence of waves of any amplitude and wavelength is sug- 


gested as an area for further study. 
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Figure 19. Plot of Roll Angle Difference Vs. Time 
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Figure 20. Plot’ of Yaw Angle Difference Vs. Time 
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Figure 21. Plot: of Yaw Rate Difference Vs. Time 
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Figure 22. Plot of Rudder Angle Difference Vs. Time 
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A listing of the FORTRAN program for simulation of 
the XR-3 is included as Appendix D. A sample input data 


deck is also included. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


This study has shown that the nonlinear mathematical 
model with six degrees of freedom adequately describes the 
motions of the XR-3 in calm water conditions. The variables 
used in the validation included thrust, longitudinal velo- 
city, pitch angle, plenum pressure, draft and vertical 
acceleration in straight runs. For validation of turning 
dynamics the variables included thrust, longitudinal velo- 
city, pitch angle, roll angle, yaw angle, yaw rate and 
rudder angle. Differences that were found to exist between 
the computer model and the testcraft were less than 5% for 
the majority of the variables compared and those differ- 
ences that exceeded 5% are considered acceptable, given 
the range of the measured variables and the accuracy of 
the sensors aboard the testcraft. 

During the course of this study, it became apparent 
that the computer model exhibits a more oscillatory, or 
underdamped, response in pitch and roll motions than does 
the XR-3 testcraft. A concurrent study by LCDR. R.A. Finley, 
USN upon this aspect of the simulation with attention to 
the modelling of the bow seal and the lift fan system 
is reported in Reference 12. 

The method of data processing set forth in this study 
to produce data in adequate form for interface with the 


Mathematical model is recommended for use in future 
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validation efforts. A considerable amount of time and 
effort was spent in signal analysis and data processing 

to prepare the measured data in a form useful for compari- 
son with the mathematical model variables. In order to 
eliminate the requirement for smoothing the re-scaled 
data, it is recommended that a study be made using a lower 
sampling rate for the A-to-D conversion to see if the high 
frequency noise can be eliminated from the digitally re- 
produced data. It is also recommended that a vibration 
study be conducted aboard the XR-3 in an attempt to deter- 
mine if the source or sources of the high frequency noise 
is the result of internal ship vibration acting on the 
Measurement transducers, gyros, and accelerometers. The 
current data acquisition system installed in the testcraft 
should be adequate for such a study. 

A very fertile area for future study would be valida- 
tion of the computer model in the presence of waves. 
Instrumentation exists to measure wave conditions at the 
XR-3 test site during measurement runs and the simulated 
waves could be input to the mathematical model. Direct 
comparison of craft behavior and simulation model calcula- 
tions could then be obtained to produce an evaluation of 


the accuracy of the simulation in sea states. 


fe 


- ors eait to tavoms sidarebs 


; er er. ioe oh a in 

ae oD < ord sisah bie 7 sed a ot 7 ne ate ¥ as 
-lisques rol Ivitess miot @ HE ats | bexvawem st req ong 

| + ~sb76 nl (eeldsitay Isdom tao tt anata se a 
s-et alt yokdtoome rt Fuenmeatypey ‘ort ots 
leu sham ad yhuta & tadd Debnowneses Slee 


id adt tl ssa af nofetevros G-of-A sit a0T ofat eee 


ij a ra +, bs : sa wee  E 
x Ys F tt BDafsenimifs $0 mas eeior 
si - ae : a Po i e ve Yr + » P 
i ¥ & + isl B Dail Wwogeat UE fs ot + + th BD 5 
‘oa. - 
} ig se I a A SLBOGE DE 5M ENOS ge 
] 
_ “<a - 7 i 
pet nto E STwvoOeB iO: Soro 
adit fro arc . ¥ gine sot to Fis 
are x 


" @i1oTe { teat 2 oliiatipos ef 
o4, ee 

~ - - * Reta 

5 raupoOneS 


= 
re 
4, 
; 
’ 


t-rot sets olivvet yas 

aly 
vod ra¥ew Bo Sorte j att ot Lebo teteqmbo aie 
r 5 a mf Bert iho 9 yivancm ot étekxe nodtepas 


» oe . 
teieente odd” by ut Greervesst aaniavt ev is oe 


| see it noite fn ony tw ‘ert 
baa 


- ul . 


ti qu je 


APPENDIX A 


LISTING OF FORTRAN PROGRAM FOR ANALOG TO 


DATA 


DIGITAL CONVERSION OF MEASURMENT 
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